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ABSTRACT 

We present 11.2 yum observations of the gravitationally lensed, radio-loud 
Zs = 2.64 quasar MG0414+0534, obtained using the Michelle camera on Gemini 
North. We find a fiux ratio anomaly of A2/A1 = 0.93±0.02 for the quasar images 
Al and A2. When combined with the 11.7 fim measurements from Minezaki 
et al. (2009), the A2/A1 fiux ratio is nearly 5a from the expected ratio for a 
model based on the two visible lens galaxies. The mid-IR fiux ratio anomaly 
can be explained by a satellite (substructure), 0'.'3 Northeast of image A2, as 
can the detailed VLBI structures of the jet produced by the quasar. When 
we combine the mid-IR fiux ratios with high-resolution VLBI measurements, 
we find a best-fit mass of io^ '^='='^-^ Mq inside the Einstein radius for a satellite 
substructure modeled as a singular isothermal sphere at the redshift of the main 
lens [zi = 0.96). We are unable to set an interesting limit on the mass to light 
ratio due to its proximity to the quasar image A2. While the observations used 
here were technically difficult, surveys of fiux anomalies in gravitational lenses 
with the James Webb Space Telescope will be simple, fast, and should well 
constrain the abundance of substructure in dark matter haloes. 



1. Introduction 



The missing satellite problem, where cosmological simulations of cold dark matter 
(CDM) predict a significantly larger fraction of lower mass satellites aroun d galaxies tha n 



is de t ected, represents a major puzzle in the study of structure formation (jKlypin et al. 



1999 



Moore et al. 



19991 ). In the case of the Milky Way, the Sloan Digital Sky Survey has 



steadily found additional, faint satellites , but the total numbers are still far lower t 
expected abundances of subhaloes (e.g.. 



lan the 



Willman et al 



2005 



Belokurov et al. 



20071 ). The 
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simplest solution in the context of CDM models is to suppress star formation in low-mass 
satellites, probably through heating and b aryonic mass loss as the universe re-ion izes to 
leave a population of dark satellites (e.g., 



Klypin et al 



1999; 



Bullock et al. 



20001 ) with 



assistance from the relative st reaming of baryons and dark matter after recombination 



( iTseliakhovich fc Hirata 



2010l ). Locally there is some hope of finding these dark satellites 



through 7-rays emitted by dark matter annihil ation, but the 



heavily on the properties of dark matter (e.g.. 



Strigari et al. 



ikeli^ 



lood of detection depends 



20081 ). 



Gravitational lensing o ffers an sl. 
lenses (e.g. MG0 414 



Kochanek et al. 



2006 



0534, 



Ros et al. 



ternate method to detect subha^ 
MG2016+112, 



2000 



Fadely fc Keeton 



Kochanek 



2012 



B1938+666, 



oes. In some strong 
HI 



2006 



Vegetti et al. 



]0435-1223, 



20121 ) ■ satellite 



galaxies of the primary lens can be detected astrometrically through their effects on the 
image positions. Potentially, the most massive satelhtes can also be detected through 



the ir effects on time delay 



by 



Mao fc Schneideil (119981 ) and 



ratio s (iKeeton fc Moustakas 



20091 ). But as emphasized 



Metcalf &: Zhaol (120021 ) . image fluxes are sensitive to 



perturbations even from very low-mass satellites (< lO^M©). Moreover, we observe 



many lenses with "flux ratio anomalies" wher e the relative image 



explained by simple, central lens ga 


axies ( 


see 


Conedon & Keeton 


20oJ: 


Yoo et al. 


2006a 


b; 



Evans fc Witt 



K^ratzer et al 



2003 



2011 



jrightnesses cannot be 



Kochanek &: Dalai 



2004 



While current results 



2002 



Vegetti et al. 



2OI2I), 



are consistent with the expectations of CDM ( iDalal &: Kochanek 
there are still too few lenses to precisely test the CDM model using flux ratio anomalies. 
For this technique, the challenge is that image fluxes can also be affected by extinction 
and microlensing in the primary lens at the most easily observed wavelengths, although 
the wavele ngth dependencie s of these three effe cts can help to disentangle them from each 



other (e.g.. 



Agol et al 



2009 



Munoz et al 



20111). The millilensing signal by substructures 



is best isolated at radio or rest-frame mid-infrared (mid-IR, ~4-100 fim) wavelengths, 
where both extinction and microlensing effects from the lens galaxy are negligible (stars 
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can only magnify sources smaller than ~0.1 pc). The mid-IR images are also unaffected 
by interstellar scattering, or scintillation, which can further perturb the radio images, 
although few lenses show the strong r adio wavelength-dependent signatures expected from 



these effects ( Kochanek fc Dalai 



120041 ) ■ The biggest problem for using radio wavelengths is 
simply that most quasars are radio-quiet, leaving only a handful of objects to be studied. 
While mid-IR observations are challenging, they are ideal for searching for substructures in 
lens galaxies. 

Here we present mid-IR observations of the quadruply lensed, radio-loud = 2.64 



known to be affected 


3y extinction and micro 


1992 




Lawrence et al. 


I995I: 


Falco et al. 


1997 



Bate et al 



2008 



2011 



near-IR ( 


Hewitt e1 


t al. 


Poolev et al. 


2012 





but the flux ratio anomaly persists into the mid-IR and radio . Our new measurements have 



slightly higher precision than those by 



Minezaki et al 



(|2009[) and we then consider more 



MacLeod et al. 



fl2009[ ). where a 



detailed models. This work is an extension of the study in 
companion lens galaxy was detected through the mid-IR flux ratios of the four-image lens 
H1413+117. Here, we show that the flux ratio anomaly in MG0414+0534 indicates the 
presence of low-mass substructure in the lens based on models of the individual or combined 
radio and mid-IR observations of this lens. Our data are described in Section |2l our lens 
models are presented in Section [3l and our conclusions are in Section |H 



Flux measurements and errors 



20041 ) on Gemini North 



We observed MG0414-h0534 with the Michelle camera (IRochd 
at 11.2 /im (F112B21 filter) on 2 Jan 2006. The observing time was 1 hour, of which 940.8 
seconds were spent on source. The data were processed with the standard Gemini pipeline 



mireduce. Our analysis starts with the coadded chop and nod subtracted image. These 
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initial images have 16 vertical stripes due to the 16 readout channels. We corrected for 
the stripes by subtracting the median of each stripe after masking the region containing 
the quasar images. This procedure also removes any residual sky flux. Figure [T] shows the 
resulting image. 

We fit the four lensed images using a 2-D Gaussian to model the PSF, fixing the 



relative i mage positions to those measured from the Hubble Space Telescope (HST) 



image4_| (iFalco et al. 



19971 ). Since the sky is so much brighter than the quasar images, 
we use the standard deviation of the background pixels as an estimate of the errors in 
the fiux of each pixel. The model has nine free parameters: the fiuxes of the four images 
(4), the position of image B (2), the standard deviation of the Gaussian PSF along the 
two principal directions (2), and the angle of the major axis. We found the best-fit model 
using non-linear Levenberg-Marquardt optimization. We used a Markov chain Monte Carlo 
(MCMC) with 10^ steps to compute the parameter uncertainties, and ran multiple chains to 
check the convergence of the best-fit parameters and their uncertainties. We also estimated 
the errors by simply computing the goodness of fit as a function of the fiux of each image, 
marginalizing over all other variables. 

Figure [1] shows the best-fit model and its residuals, and Figure [2] shows the MCMC 
and Ax^ estimate for the fiuxes (in arbitrary units) and their uncertainties. In order to 
include any covariances between the image fiuxes, we use the actual distribution of fiux 
ratios from the MCMC approach to estimate the fiux ratio uncertainties. We calibrated 
the fiux measurements using observations of the standard star Rho Orionis (HR 1698) and 



the 



Cohen et al. 



( 1l999l ) fiux calibration at 11.2 fim (the interpolated fiux is ~ 6.43 Jy). 



We obtained a total fiux of 33.4 ± 1.3 (statistical) ±3.2(systematic) mJy at 11.2 nm frorn 



aperture photometry, which is lower than the value of 39.2 ± 1.4 mJy in 



Minezaki et al 



-"^http: / /www. cfa.harvard.edu/castles/ 
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Fig. 1. — The left, middle, and right panels show the 11.2 /im data, the best-fit model based on 2-D 
Gaussians, and the model residuals, respectively, for MG0414+0534. In the bottom panels, the images have 
been scaled to increase the visibility of B and C. 



(120091 ). This discrepancy may be due to the fact that our filter has a shorter central 



wavelength and that the SED is rising towards longer wavelengths. The resulting mi d-IR 



flux ratios are listed in Table [T] along with the mid-IR results from 



Minezaki et al 



(120091) 



and the combined results. The results agree given their mutual uncertai nties and show 



20081 ). persists into the 



that the anomaly, while not as strong as in the near-IR (iBate et al.l 
mid-IR. This wavelength behavior is a first indication that the anomalous flux r atio in the 
near-IR is not produced solely by microlensing, as assumed in 
by substructure in the lens. 



Bate et al.i (i2008l) . but also 
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Fig. 2. — MCMC probability distributions for the MG0414+0534 image fluxes (in arbitrary units). The 
vertical dashed lines show the mean and la error bars determined from the best-fit Gaussian curve. These 
are compared to the best-fit flux value and 68% confidence limits (or parameter values at which Ax^ — 1 
from the minimum x^) that result from the minimization technique (see text), which are shown by vertical 
dotted lines. The Gaussian widths (a) are listed in the top right corners. For the uncertainties reported in 
Table [H we instead use the dispersion of flux ratios from the MCMC. 



Table 1. Relative Mid-IR Fluxes for the Images of MG0414+0534 



Images 




Flux Ratio 






This work 


Minezaki et al. (2009) 


Combined 


A2/A1 


0.926 ± 0.025 


0.90 ± 0.04 


0.919 ± 0.021 


B/Al 


0.338 ± 0.017 


0.36 ± 0.02 


0.347 ± 0.013 


C/Al 


0.145 ± 0.016 


0.12 ± 0.03 


0.139 ± 0.014 
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Models 



In a multiple-image strong gravitational lens, a close pair of bright images such as 



Al and A2 results if the quasar is near a fold caustic created by the main lens ga. 



Keeton et al 



axy. 



20051), 



Their flux ratio is expected to be unity for a smooth lens model (e.g 
whereas we measure a flux ratio of A2/A1 = 0.93 ± 0.02 for MG0414+0534. Since we are 
insensitive to extinction or microlensing effects created by the lens galaxy in the mid-IR, 
the anomalous ratio must be attribu ted to millilensing by substructure in the lens, as was 



also argued by 



Minezaki et al. 



(120091 ) 



.1 



200l[ ). using the combined mid-IR 



We model the system with LENSMODEIJfl flKeetonI 
flux ratios listed in Table [T] (last column) as constraints on our models. We use the measured 
HST WFPC/WFPC2 po s itions of the ma in, central lens galaxy, the secondary lens "X" 



( ISchechter fc Moore 



database (Falco et al. 



993 



Ros et al. 



20001 ). and the image positions from the CASTLES 
19971 ). The component positions are listed in Table |21 The main lens 
galaxy (Gl) is modeled as a singular isothermal ellipsoid (SIE) combined with an external 
shear, as this is the simplest plausible model for a lens. There is considerable evidence in 
favor of the isothermal profile, and the properties of 4-i mage systeins other than time delay s 



Kochanek 



2006 



Koopmans et al. 



are v ery insensitive to the radial mass distribution (see 
20091 ). We include weak priors on the ellipticity of the main lens (1 — 6/a = ± 1), and the 
external shear (7 = 0.05 ± 0.05). The expectation value of for the ellipticity is based ori 



the r elatively round isophotes derived from the HST photometry of the lens (IFalco et al. 



19971 ). For object X (hereafter "G2"), we tried two different models: a singular isothermal 



sphere (SIS) or a pseudo-Jaffe (pJaffe) profile with a scale length a' that is allowed to vary. 



Using the mid-IR flux ratios as constraints, we obtain the best-fit two-galaxy models 



^http: / /redfive. rutgers.edu/ ~keeton/gravlens/ 
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Table 2. MG0414+0534 Image and Lens Positions 



MG0414 component HST position (") Rotated 0.1° (") Error (") 



Al (-0.600,-1.942) (-0.596,-1.943) 0.003 

A2 (-0.732,-1.549) (-0.729,-1.550) 0.003 

B (0.000,0.000) (0.000,0.000) 0.003 

C (1.342,-1.650) (1.345,-1.648) 0.003 

Gl (0.472,-1.277) (0.474,-1.276) 0.003 

02 (=X) (0.857,0.180) (0.857,0.181) 0.011 



Note. 



(IFalco et al 



The positions are taken from the CASTLES database 



19971 ) and are relative to image B, where negative R.A. values 
are eastward of image B. We rotate the positions counter-clockwise by 0.1° 
when adding the radio constraints to our lens models. 
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described in Table |3l The models have 17 total constraints (the positions of 4 lensed images 
and both lens galaxies, the 3 flux ratios, plus the priors on the ellipticity e and shear 7) 
and 12 free parameters (the mass scale bi, position, ellipticity e and position angle 6^ for 
Gl, the amplitude 7 and position angle 6^ of the external shear, the mass scale 62 and 
position for G2, and the source position), plus one more free parameter for the scale length 
a'g when using a pseudo-Jaffe profile. In the Table, we describe our best two-galaxy models 
when using SIE+SIS or SIE+pJaffe profiles. We find a slightly better fit when using an 
SIS profile for G2, with = 31.1 for dof=5 degrees of freedom and a best-fit fiux ratio 
A2/A1 = 1.01. The discrepancy between the observed and predicted mid-IR fiux ratios 
in our two-galaxy models suggests the presence of substructure along the line of sight or 
associated with the main lens galaxy at nearly a 5a level. 

We modeled the substructure by adding a third lens galaxy (G3) with an SIS profile 
at various positions surrounding the lens and searched for a better fit to the observed 
fiux ratios. Here, we adopt an SIS rather than a pseudo-Jaffe profile for G2 since the 
former yielded a better fit among the two-galaxy models. At each position in the grid, 
the mass scales bi, 62, and 63 for galaxies Gl, G2, and G3 were allowed to vary, as well 
as the ellipticity of Gl and the external shear, until a best-fit three-galaxy model was 
obtained. Here, dof=6, since we are holding the positions for Gl and G2 fixed for each 
model. Figure [3] (left panel) shows the resulting Ax^ (relative to the two-galaxy model). 
It can be seen that the addition of substructure improves the two-galaxy model when 
placed Northeast or Northwest of image A2, although the range of positions that lead to a 
similar improvement is more limited Northwest of A2. We use this approximate location 
for G3 (Northeast of A2) as a starting point in our three-galaxy model and obtain the 
best-fit parameters in Table IHl where the coordinates of G3 are treated as free parameters. 
We find that an SIS with 63 ^ 0'.'0028 located at {x^^ys) = (-0'.'91, -1'.'56) leads to the 
best improvement, with Ax^ = 25 (x^/dof = 3). Assuming the measured redshift of 
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Table 3. MG0414+0534 Modeling Results 





Gl+G2(pJaffe) 


G1+G2 


G1+G2+G3 


G1+G2 


G1+G2+G3 




(Mid-IR only) 


(Mid-IR + radio) 


X /QOt 


28.8/4 


31.1/5 


6.11/2 


103/41 


ov.o / 00 


^2 
Apos 


4.83 


5.46 


0.142 


8.45 


5.51 


v2 


22 7 


9'^ Q 


A 71 


on ^ 
yu.u 


52 1 


v2 
AgaZ 


0.746 


1.18 


0.539 


3.17 


1.06 


^2 


0.613 


0.634 


0.720 


0.928 


0.849 


01 


1'.'104 


1'.'068 


1'.'091 


1'.'083 


1^'084 




0.298 


0.302 


0.224 


0.242 


0.238 


f) C°R of N^l 


-88.5 


-88.0 


-85.4 


-84.5 


-83.9 


1 


0.086 


0.087 


0.091 


0.097 


0.094 


e-y (°E of N) 


47.5 


47.2 


54.7 


53.6 


53.1 


62 


0'.'199 


0'.'193 


0'.'170 


0'.'185 


0^'176 




2'.'53 










63 






0'.'0028 




0^007 ± O'.'OOl 








-0^'91 




-0'.'97 


ys 






-1^56 




-1'.'39 


Mid-IR Flux Ratios 












A2/A1 


1.0018 


1.0058 


0.9218 


1.0320 


0.9227 


B/Al 


0.3318 


0.3415 


0.3438 


0.3308 


0.3354 


C/Al 


0.1690 


0.1712 


0.1691 


0.1614 


0.1680 



Note. — The best-fit parameters for our two- and three-galaxy models based on 
the mid-IR data alone (left three columns) and after adding the VLBI constraints 
(right two columns). In the first column, G2 is modeled as a pseudo-Jaffe potential 
with scale length a'2. In the remaining columns, G2 is modeled as an SIS. The 
subscripts differentiate between the main lens (1), secondary lens galaxy G2/object 
"X" (2), and a third lens mass (3). The is broken down into the contributions 
from the image positions (Xpos); ^^^x ratios (X/zux)' ^^^^ position (Xgai)' 
priors on e and 7 (xLior)- 
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zi = 0.96 for Gl (iTonry &: Kochaneklll999l ). we estimate a best-fit mass of ~ 3 x 10^ Mq 



and ~ 10^*^ Mq enclosed within the Einstein radii of G3 and G2, respectively, based on the 
mid-IR constraints alone, as compared to ~ 5 x 10^^ Mq for Gl. The statistical significance 
of finding Ax^ = 25 for 3 additional parameters is 0.28 (a 72% confidence level) based on 
the F-test. When G2 is modeled as a pseudo-Jaffe potential, we find similar results for a 
three-galaxy model but with x^/dof = 5.8 (and an F-test confidence level of 45%). We 
conclude that G2 is better modeled as an SIS than a pseudo-Jaffe potential and therefore 
do not list the model parameters for a three-galaxy (SIE+pJaffe+SIS) model here. 
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Fig. 3. — Ax^ improvement resulting from adding a third lens galaxy to MG0414+0534 as a function of its 
position. In the left panel, only the mid-IR constraints are used in the models. In the right panel, only the 
radio constraints are used. The white regions are where the either matched or exceeded the best value 
from the 2-galaxy models. North is up, East is left. 
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3.1. Radio Constraints: 8.4 GHz VLBI Structures 



Fortunately, we can further test this model 



by adding the flux ratios and positio ns of the 



2000 



Trotter et al 



20001 ) ■ where 



components detected in the radio with VLBI ( iRos et al 
each quasar image is resolved into core and jet structures. We added the core components a 
and b detected in the VLBI maps, and the VLBI jet structures c and d . The radio position s 



Trotter et al. 



fcoooh 



were constrained using the error ellipses shown in Figure |H Note that 
constrained their models using tighter constraints on the component positions while ours 
are more conservative, especially along the jet structures, to allow more freedom in the 
component cross-identifications. We find a better agreement (Ax^ = 3) between the optical 
and radio core positions when shifting the HST images counter-clockwise by 0.1°, which is 
well within the HST astrometric uncertainties. For the remaining models, we assume the 
rotated HST positions for the mid-IR images (gray circles) and lens galaxies. Our best-fit 
component positions are indicated with symbols. 

We list the best-fit model parameters for a two-galaxy model in Table [31 Here, we 
adopt an SIS profile for G2 since it provides a marginally better fit in this case than a 
pseudo-Jaffe profile. Our two-galaxy model has dof = 41, since there are four additional 
sources, each with three fiux ratios and four pairs of (a;, y) image coordinates as constraints, 
and four additional source positions as free parameters. 

Figure [5] shows the resulting improvement when adding a third lens galaxy in the 
same fashion as before, but including the constraints from all four radio sources (a similar 
map using the radio constraints only is shown in the right panel of Figure [3]). The best-fit 
mass scale 63 is also shown. The overall results are consistent with those based only on 
the mid-IR data, but with a narrower region extending to the Northeast (145° North of 
West) for the position of the third galaxy. The parameters for the best three-galaxy model 
with (x3,|/3) = (— 0'.'97, — 1'.'39) and x^/dof = 1.6 are hsted in Table [3l If we vary the mass 
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parameter of G3, we find it is well-constrained to 63 = 0'/007 ± O'.'OOl despite the freedom 
to simultaneously adjust all other model parameters, including the external shear. This 
corresponds to IQ'^-^^^-'^ Mq enclosed within the Einstein radius at zi = 0.96. The statistical 
significance of this result is > 99% based on the F-test. 

The critical and caustic curves for our final three-galaxy model are shown in Figure El 
The third galaxy distorts the critical curves near images Al and A2, bringing their fiux 
ratio back into agreement with the mid-IR data. G3 may be associated with a faint dwarf 
galaxy or a nonluminous darkmatter substructure, as there are no nearby objects detected 
in the H-band image which are coincident with its predicted positions. Unfortunately, in 
the H-band image where we would best be able to detect G3, the estimated position lies 
almost exactly on the first Airy ring of image A2 (where the image subtraction residuals 
are high) and in a region with substantial emission from the quasar host galaxy, making it 
impossible to place useful limits on the mass-to-light ratio of G3. 

The orientation of the external shear in our three-galaxy model, indicated in Figure Ej 



is in the same general direction as ano ther object 



the H-band image. As pointed out by 



Falco et al. 



(G61 



bund 4'.'4 Southwest of image C in 



(Il997l ). G6 is a possible compact group 



of galaxies, as it consists of three visible galaxies within 1" of each other. If we assume that 
G6 and G2 are at the same redshift as Gl and scale the Einstein radii of these objects by 
their I-band fiuxes assuming a standard Faber- Jackson relation {L oc a^), we can estimate 
the defiection scale of G6 as = bG2{LQ/ Lg2Y^'^, using the fact that b oc cr^ for a n SIS . 



Falco et al 



fll997h . 



Adopting the I-band magnitudes of G2 (24.769) and G6 (23.16) from 
we estimate fog = 0'.'34. The shear caused by G6 is then ■jq = bQ/2rQ = 0.03, where rg 
is the angular separation from Gl, and the am plitude of any higher order perturbations 



is b^bci/rl = 0.013 (see iKochanek et al.l 120061 ). Since 75 < 0.09, G6 does not appear 



massive enough to fully explain the external shear in our three-galaxy model, although it 
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likely contributes some shear. The inclusion of G6 as a fourth lens mass yields a best-fit 
with be = (y.'21, a reduction in shear amplitude (0.061), and a minor change in the shear 
orientation (53.3°) with xVdof = 52.6/37 = 1.4. 

In an attempt to bring our reduced closer to unity, we tried adding a fifth lens 
galaxy on a grid of positions while optimizing the overall lens properties. We find the best 
improvement (Ax^ ^ 7) near image C, extending to a region ~2" North of G6. However, 
the reduced does not improve (46/34 = 1.4). We conclude that additional substructures 
are unlikely to lead to a significant improvement of the fit given the constraints adopted 
here. 
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Fig. 4. — The constraints on the positions of the five sources used in our final model are shown 
in the image plane as error ellipses . Our esti r nates for the VLBI core components (a and b) 



and jet structures (c and d) from 



Rosetal 



( I2OOOI ) are shown as black ellipses. The error 



circles for the mid-IR positions are shown in gray. The best-fit positions for each component 
are indicated by the symbols and are consistent with each error ellipse. The larger symbols 
correspond to the mid-IR image positions, and the smaller, thick gray symbols correspond 
to VLBI component b. 
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Fig. 5. — The left panels show the A^^ improvement resulting from adding a third lens galaxy as a function 
of its position, as in Fig.[3l but using both the mid-IR and radio constraints on the models. The right panels 
show the best-fit critical radius for the substructure. The bottom panels give an expanded view of the region 
indicated by the box in the top-left panel. The positions of two nearby galaxies G4 and G5 seen in the 
H-band image are also indicated. 
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Fig. 6. — Best-fit critical (thin) and caustic (red) curves for our tfiree-galaxy model when 
including both the mid-IR and radio constraints. Galaxy (image) positions are indicated 
with open circles (crosses). The direction of the external shear is indicated with arrows, and 
points toward another galaxy (G6) detected in the H-band image. 
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4. Discussion 



Gravitationally lensed quasars are useful astrophysical tools because they allow us to 
study both the lens mass distribution and the background quasar. Unlike microlensing, 
where time variability can be used to isolate the effect, the effects of millilensing need to 
be disentangled from other physics that modify flux ratios. Mid-IR observations are ideal 
for this purpose becau se the effects of both microlensing and extinction are negligible. In 



MacLeod et al. 



(120091 ). we used the mid-IR flux ratios for the lens H1413+117 to show 
that flux ratio anomalies can detect companion lens galaxies, albeit a luminous companion 
in that system. This technique can also be used for detecting d ark subhaloes in the lens. 



Minezaki et al. 



down to lO^Mr;. flMao fc Schneider 



1998 



Metcalf fc Zhao 



20021 ). As shown here and in 



(120091 ). the quadruply lensed quasar MG0414+0534 exhibits a flux ratio 



anomaly at mid-IR wavelengths that suggests there is millilensing substructure in the lens 
or along the line of sight to the quasar. 

When we model either the mid-IR or radio data or both using only two lens galaxies 
(the primary lens Gl and a secondary lens G2/object X), we are unable to reproduce the 
observed flux ratios. When we add a third lens G3 near image A2, we obtain a signiflcantly 
improved flt to the flux ratios. The optimal location of the substructure, roughly 0'.'29 to 
the Northeast of image A2 and 1'.'45 East of Gl, was determined by adding a third lens 
on a grid of positions in the lens plane and varying its mass along with the overall lens 
properties. We are unable to associate the substructure with any source in the H-band 
image from HST, which suggests the third galaxy is either very faint or nonluminous. 
We derive a mass of iO''-^='=°-^ Mq enclosed within the Einstein radius at zi = 0.96 for 
the substructure when modeled as an SIS, and our results are statistically signiflcant at 
> 99% confldence based on the F test. Our flndings are generally consistent with the less 



detailed models of 



Minezaki et al. 



(120091 ). The anomaly in the mid-IR suggests the presence 
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of a third lens galaxy which we assume is at the redshift of the main lens, alt hough the 



struc t ure could theoretically li e at some other redshift along the line of sight fIXu et al. 



2012 



Inoue fc Takahashi 



20121 ). There is also evidence for some contribution to the flux 



ratios from the neighboring galaxy group (G6) in the H-band image. 

Besides being relevant to studying the background quasar, the search for subhaloes in 
the lens is also cosmologically important, because it offers the only means of identifying 
dark, low mass halos other than detection of 7-ray annihilation signals by Fermi. Sufficiently 
massive structures can be detected throu gh their perturb ations of image po sitions as seen 



with object X here, or in MG2016+112 (IKochanek 



20061 ) and B1938+666 fIVegetti et al. 



2OI2I ). Lower mass structures can only be detected through flux ratio anomalies, where 
the challenge is to obtain flux ratios unaffected by absorption or microlensing. The mid-IR 
satisfies these requirements but is limited at present by either poor resolution, as in Spitzer 
observations, or poor se nsitivity, as in gr ound-based measurements. The LBTI mid-IR 



imaging interferometer ( iHinz et al 



20031 ) on the Large Binocular Telescope (LBT) is a 



promising near-term solution, and a complete survey will be trivial with the advent of the 
James Webb Space Telescope. 



This work is based on observations obtained at the Gemini Observatory (with Program 
ID GN-2005B-Q-43), which is operated by the Association of Universities for Research in 
Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (United States), the Science and Technology 
Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), Ministrio da Cincia e Tecnologia 
(Brazil) and Ministerio de Ciencia, Tecnologa e Innovacin Productiva (Argentina). This 
work is based in part on observations made with the NASA/ESA Hubble Space Telescope. 
Support for program GO-7495 was provided by NASA through a grant from the Space 
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Telescope Science Institute, which is operated by AURA, Inc., under NASA contract 
NAS5-2655. EA is partially supported by National Science Foundation CAREER Grant 
No. 0645416. CSK is supported by NSF grant AST-1009756. 

Facilities: Gemini (Michelle), HST (NICMOS). 
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